sympy2c: From symbolic expressions to fast C/C++ functions and ODE solvers in Python

Computer algebra systems play an important role in science as they facilitate the development of new theoretical models. The resulting symbolic equations are often implemented a compiled programming language order to provide fast and portable codes for practical applications. We describe sympy2c, Python package designed bridge gap between numerical implementation model. sympy2c translates SymPy C/C++ code that is optimized using transformations. functions can be conveniently used extension module Python. within PyCosmo solve Einstein–Boltzmann equations, large system ODEs describing evolution linear perturbations Universe. After reviewing functionalities usage we its optimization strategies. This includes, particular, novel approach generate ODE solvers making use sparsity Jacobian matrix. demonstrate performance test case. general potentially useful various areas computational physics. publicly available at https://cosmology.ethz.ch/research/software-lab/sympy2c.html.